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In this paper we examine the relationship between the loss of the 
TEoi mode in a multimode circular waveguide and the waveguide ge- 
ometry. First-order perturbation theory solutions of the coupled line 
equations were used to predict the loss from measurements of the 
waveguide geometry. The predicted loss disagreed with the measured 
loss for the 14 km long waveguide line in the WT4 field evaluation test. 
An analysis of the coupled line equations which considers the effects 
of both first- and second-order mode conversion is described. The ap- 
propriate coupling coefficients are derived and discussed. Second-order 
perturbation theory solutions of the coupled line equations are devel- 
oped and examined. Losses predicted from the measured waveguide 
curvature and the second-order theory agree well with measured loss 
data. The results indicate that axial curvature leads to significant mode 
conversion-reconversion effects between the TEqi-TMu, TEq\-TE\2, 
TM n -TM 2h TE12-TM21 modes in the installed waveguide of the WT4 
field evaluation test. The TM 2 \ conversion is due to curvature fluctu- 
ations with a characteristic wavelength of 10 to 30 m. 

I. INTRODUCTION 

The transmission characteristics of the TE i mode in circular wave- 
guide are quite different from the characteristics of cable or wire 
transmission media currently used in the Bell System. The waveguide 
medium has the unique property that its heat loss decreases with in- 
creasing frequency. Heat losses 1 as low as 0.3 dB/km are obtained at 110 
GHz in 60 mm diameter waveguide. However, since the diameter of the 
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waveguide is much larger than the operating wavelength of 2.75 mm at 
110 GHz, several hundred spurious modes can also propagate, in addition 
to the TE i mode. The spurious modes have slightly different phase and 
group velocities than the TEoi mode. Spurious modes are excited in 
practice because the waveguide is not a geometrically perfect right cir- 
cular cylinder, and there is a continuous exchange of energy between the 
TEoi and spurious modes. This mode conversion-reconversion results 
in added loss and delay distortion. A detailed understanding of the re- 
lationship between spurious mode generation and waveguide geometry 
is required since even small distortions in the guide cross section (of the 
order of 1 /im) can cause substantial loss. 2 

The mode conversion problem can be modeled by an infinite system 
of coupled transmission line equations 3 for which approximate solutions 
have been obtained by first-order perturbation theory. 4 At the beginning 
of the WT4 field evaluation test the first-order theory was considered 
to be a valid and accurate method for predicting the TEoi mode con- 
version loss. Detailed measurements of the waveguide cross-section 
geometry were carried out in the laboratory on 9 m long individual sec- 
tions of waveguide. The curvature of the waveguide axis, in two or- 
thogonal planes, was also measured over the 14 km field evaluation test 
route. The curvature and cross-section geometry data and the first-order 
perturbation theory solutions were used to predict the TEqi mode con- 
version loss. Waveguide curvature was found to be a significant source 
of TMn and TE12 mode generation with a predicted TEoi mode con- 
version loss of approximately 0.2 dB/km at 110 GHz. All other predicted 
first-order mode conversion effects resulted in an additional TEoi loss 
of less than 0.03 dB/km at 110 GHz. 

The measured loss disagreed substantially with those first-order 
predictions. The measured loss indicated a total TEoi mode conversion 
loss of approximately 0.5 dB/km at 110 GHz. 

In this paper, we show that the generation of the TM21 mode is a sig- 
nificant source of TEoi mode conversion loss in curved waveguide. The 
TM21 mode is not coupled directly to the TEoi mode and therefore is 
ignored in the first-order perturbation theory. It is excited via the TMn 
and TE12 modes and hence can be thought of as a "second-order" phe- 
monenon. However, as we shall see, the energy coupled into the TM21 
mode exceeds the energies coupled into the "first-order modes" TMn 
and TE12 under certain conditions. Therefore the term "second order" 
is used to describe any theory which includes the TM 2 i mode, only in 
the sense that it is a mode which is indirectly excited from the TEoi 
mode. 

An approximate model and a second-order perturbation theory is 
developed for the TM21 conversion in curved waveguide, which, in 
conjunction with the first-order theory for TMn and TE12 conversion, 
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yields results which agree well with measured results. The source of the 
TM21 generation is identified as waveguide curvature with low spatial 
frequency components with a characteristic wavelength of 10 to 30 me- 
ters. 

II. PROBLEM FORMULATION 

Maxwell's equation in the deformed waveguide may be transformed 
to an infinite set 3,5 of coupled transmission line equations expressed in 
terms of the normal modes of the undeformed guide: 

— A m (z)= L K min {z)A n (z) -jh m A m (1) 

The unknowns in eq. (1) are the normal mode amplitudes A m while h m 
are the normal mode propagation constants in the undeformed guide. 
The sum over n extends over an infinite number of modes, and the 
coefficients K m>n in the infinite system of equations consist of the 
product of a geometry deformation factor and a normalized coupling 4 ' 5,6 - 7 
coefficient C mn . The derivation of the coupled line equations and the 
normalized coupling coefficients, as well as some notation for normal 
modes, is outlined in the appendix. 

We shall use single subscripts m and n when referring to a general 
normal mode. A double subscript pn is used in referring to a specific 
normal mode such as TE P „ or TM pn . Here the p refers to the order of 
the azimuthal variation (cos p<f> and sin p<f>) while the n is the order of 
the radial variation J p (x n r). The modes in dielectric-lined waveguide, 
the principal transmission medium in the field evaluation test, are not 
strictly 6 ' 7,8 "transverse electric" TE or "transverse magnetic" TM but 
we will retain the notation. The axis of the installed waveguide is a curve 
in three-dimensional space. The polarizations of the locally excited 
spurious modes are related to the orientation of the osculating plane, 
and since this is continuously changing, the polarizations of the modes 
are also changing continuously along the waveguide axis. The net spu- 
rious mode amplitude is the sum of the locally generated modes, and to 
determine this sum it is necessary to resolve the local polarization into 
two orthogonal polarizations. Furthermore, it has been shown 4 that the 
same results are obtained by using the curvature of the projections of 
the waveguide axis on two orthogonal planes. We shall use superscripts 
t and —*■ to denote the two distinct polarizations, e.g., TE p „ or TE^,. We 
shall use TE pn with no superscripts in discussions which apply for both 
polarizations. The two polarizations for the TMji and TM^ modes are 
given in Fig. 1. Note that the plane of curvature determines which po- 
larization is generated. The waveguide installed for the field evaluation 
test has independent curvature characteristics for the two planes, hor- 
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Fig. 1 — Normal mode electric fields in dielectric-lined waveguide. 

izontal and vertical, shown in Fig. 1, and the polarization of the spurious 
mode must be carefully accounted for when estimating the TEoi loss. 
Previous analyses 4 assume K mi1l {z) to be sufficiently small so that the 
total loss can be obtained by a superposition of a number of two mode 
problems in each of which the TEoi mode and one spurious mode is 
considered. As a specific example, the TEoi-TMn mode conversion loss 
due to axial curvature c(z) in a single plane may be determined from 



dzlAni' ; |_0 /inJUnJ 



+ C„ r ][£]«« *> 



An approximate 4 expression for the expected TEoi-TMn mode 
conversion loss (all) in eq. (2) is 



(* 1 o\ = -C 2 01M S c (APoi,ii/2*) 



(3) 



Here S c is the spectral density function for c and the differential prop- 
agation constant Aj8bi.ii is defined as 

A/3oi f n = h i - fin (4) 

Expression (3) is valid in loss-free waveguide, real h m , and for long 
lengths of line with a slowly varying power spectrum, S c . 4 It is an exact 
solution for (2) in the case of a white power spectrum. 9 Computed power 
spectral density functions of the installed WT4 field evaluation test 
waveguide 2,10 diameter, curvature, ellipticity, and higher-order geometric 
distortions were substituted into the equivalent of (3). The resulting 
values of loss indicated that approximately 90 percent of the total TEoi 
mode conversion loss, aoi C , was accounted for by TMji, TMn, TE{ 2 , TE^ 
conversions due to the curvatures c# and cy of the waveguide axis in the 
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horizontal and vertical planes, respectively. Thus the predicted mode 
conversion loss was expected to be 

«oT = \ Cfcii [S cv (A^oi,ii/27r) + S C|I (A|8oi,ii/2t)1 

+ \ CS1.12 [Sc V (Aj8bi,i2/2x) + S ch (A^oi,12/2t)] (5) 

The mode conversion loss predicted by (5) differed significantly from 
the measured loss 10 as discussed in greater detail in Section IV. Moving 
piston measurements 11 in regions of the field evaluation test line with 
high curvature indicated the presence of significant levels for the TM21 
mode as well as TMn and TE12 modes. 

The measured results implied that second-order conversion processes 
of the form TE01 ** TMn ** TM21 are also significant and must be 
considered along with the first-order processes, TE01 ** TMn, in esti- 
mating the total TE01 loss. Instead of considering the two-mode model 
in (2), we must examine a three-mode model of the form 
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In (6) we are considering only one of the possible TM21 conversion paths, 
as the curvature c(z) is assumed to be in a single plane. The other paths 
are given in Fig. 2. The coupling model in Fig. 2 leads to the following 
set of coupled transmission line equations for the determination of TE01 
mode conversion loss due to curvature of the waveguide axis: 



- [A] = -j[h][A] + jc H (z)[C H ][A] + jcv(z)[C v )[A] (7) 

dz 



where 



A = [A 01 ,A\i,A\2,ALAn,An,An] T (8) 

Here the A m 's are the transmission line voltages corresponding to the 
TE01, TMji, . . . , TM^ modes. The propagation constants hoi, Mi> • • • » 
K£\ of the modes in straight waveguide are the elements of the diagonal 
matrix [h]. Note that /i^, = h^ and that no superscript is required for 
the TE01 mode. Corresponding differential propagation constants, A# m>n 
are given in Fig. 3. 
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Fig. 2 — Curvature coupling mechanisms. 



The coupling coefficient matrices [Ch] and [Cy] consist of the nor- 
malized curvature coupling coefficients C m , n from mode m to mode n. 
They may be determined from eq. (57) in the appendix. The normalized 
curvature coupling coefficients in 60 mm diameter guide, with a 180-/im 
polyethylene liner are shown in Fig. 4. [Ch] and [Cy] must be multiplied 
by the local horizontal or vertical curvature, c//(z) or cy(z), respectively, 
to obtain the local coupling coefficients. From Fig. 2 we have 
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Fig. 3— Differential propagation constants in 60-mm lined waveguide. 

Prediction of TE i loss based on the set of coupled line equations in (7) 
agrees well with measured results as discussed in Section IV. 

III. SOLUTIONS— COUPLED LINE EQUATIONS 

The set of seven coupled line equations in (7) has been integrated 
numerically for sets of measured curvature data up to 800 m in length. 

TE i TRANSMISSION WITH AXIAL CURVATURE 1855 



35 



E 20 



UNITS 


K 0I .2,/10 


C mn (DIMENSIONLESS) 


/ 

1 
1 




f 01.21 ' m ' 


1 
1 
1 






1 s 


£, = 2.28 


1/ 


t = 180/J.m yi 


r c oi.i2 


/I 




x I 




/ 1 




/ 1 




/ 1 






Coi.ii 


/ i 




/ i 






rGyigx 


/ J/s' 




/ ^^^^ 




/ /// 




/ yS\r / 




/ S^r * 




/ s's' t 




/ s^S^ t 




/ X^s^ t 




/ >>^ / 




/ s£-^ / 




/ y^^ ' 




S ^^ f 




/ ^^ ' 




^r / 




^"/ ' 




*\s ' 




X / 




r / 




s 




y 


C 12.21 


■ y~ 


















1 1 1 1 





40 



80 100 

FREQUENCY, GHz 



120 



Fig. 4 — Coupling coefficients for 60-mm-diameter lined waveguide. 



The numerical method yields solutions with an accuracy of better than 
1 percent. The method consists of a direct numerical evaluation of the 
transmission matrix for incremental sections 5 cm long. Direct multi- 
plication of the individual matrices yields the transmission matrix for 
the desired sections of line. The results agree well with measured data 
as discussed in Section IV and are a confirmation of the validity of the 
mode conversion model in Fig. 2. The numerical integration must be 
carried out at a number of discrete frequencies to determine the loss 
characteristics over a given frequency band and is an expensive process, 
since 10 minutes or more of computer time may be required for a 500 m 
length. Thus approximate solutions are required. 
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For the case of weak coupling, the seven-mode system of equations 
in (7) can be reduced to a superposition of three-mode equations similar 
to (6). It is convenient to define a normalized mode amplitude G m in 
terms of A m 

Am = e-#**Q m (11) 

On applying Picard's method 4 of successive approximations to (6), we 
find that On and G21 are given to first- and second-order respectively 
by 

Gil = J fi C i,iic(s)e-JMoi,ns ds (12 ) 

G21 = " f * Cn a cWe-^«" f S C01 iic(t)e-i*tow'dt ds (13) 
Jo Jo 

For slowly varying curvature, c{z), it can be shown that 

Q*~j £ \c n ,2ic(s) (^jj^c(s)}\e-Mw»d8 (14) 

There is a simple physical interpretation for the terms in (14) which 
we will examine instead of presenting the details of the analysis. As the 
TE01 mode travels through a region of waveguide with slowly varying 
curvature in a single plane, the local field structure for the TE01 mode 
is slightly distorted from the structure in the straight waveguide. The 
degree of distortion can be quantified by expanding the fields in the 
curved guide in terms of the normal modes for the straight waveguide. 12 
In doing so, it is found that the TMn component has a magnitude given 
by 



Woi,ii / 



which is one of the terms in (14). The TMn component travels with the 
same phase velocity as the TE01 mode and couples to the TM21 mode 
with a coupling of strength Cn,2i c(z). This is the second coupling term 
in (14). 

Thus the net coupling between the TE01 and TM21 modes in wave- 
guide with slowly varying curvature, c, is 

Coi,iiCu2i 9 , . 

— tl c 2 (z) 

Aj8oi.il 

Equation (14) is similar in form to (12) and thus the three-mode 
mechanism, TE01 ** TMn ** TM21, has been converted to an equivalent 
simple two-mode mechanism TE01 ** TM21 with an effective normalized 
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coupling coefficient, C i,2i, given by 

r C01.11C1U1 n -v 

^01,21 - KQ \ LO > 

Aj8oi,n 

The waveguide geometry enters (14) via the factor c 2 (z). The reader is 
cautioned against confusing (15) with the normalized coupling coefficient 
between the TE i and TM21 modes due to ellipticity of the wave- 
guide. 

The formula for the expected value of the TE01 ** TM21 mode con- 
version loss, a$i, is similar to (3) and is given by 

«& = ^ lt21 S c2 (Aft>i,2i/27r) (16) 

The loss in (16) is proportional to the spectral density function for the 
square of the waveguide curvature at a spatial frequency of A0oi,2i/2ir 
c/m. From Fig. 3 it may be observed that A/3oi,2i/2?r is less than 0.1 c/m 
for frequencies greater than 80 GHz, and thus curvatures with long 
wavelengths, 10m or more are the source of TE01-TM21 mode conversion 
loss in (16). 

In Ref. (10) it is shown that the curvature power spectral density for 
the field evaluation test waveguide in the to 0.1 c/m region of the 
spectrum is three orders of magnitude larger than the spectrum in the 
0.3-1 c/m region which causes TMn, TE12 conversion. Thus the as- 
sumptions in previous analyses 4 of a slowly varying spectrum are clearly 
not true and the existence of additional mode conversion loss compo- 
nents is not surprising. 

In general, for curvature in two planes, we must consider all of the 
mechanisms in Fig. 2 to determine the total TE01 loss. It is convenient 
to define an equivalent normalized coupling coefficient, if 01,21 

tfoui s 2 [ c °MiC",tt + S«ag*a ] (17) 

L A/Jn A0i2 J 

where 

Zhx = Afa ' 2 ' + Ag "' ' , Afe - Afe ' 21 + Afe - 01 (18) 

z z 

A plot of Koi,2i vs. frequency is given in Fig. 4. 

We then find that the mode conversion due to curvature may be ap- 
proximately treated as a superposition of the two-mode process shown 
in Table I. 
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Table I — Mode conversion due to curvature 

Normalized 
Conversion Geometrical coupling 

process factor coefficient A/J Loss 

TEoi-TMl, c H Coi.il A0bi.ii 5 CSluSc!, (A0bi.ii/2x) 

TEoi-TMH c v Coi.il Aj9bi.ii \c 2 olM S cv (A0 O i.n/2r) 

TEoi-TEl 2 c h C01.12 A801.12 sC&, M S eJ ,(A/3bi.i2/2ir) 

TE01-TE15 cv C01.12 A0O1.12 g C oi.i2Sc V (A0oi,i2/2ir) 

TEoi-TM 21 cft-c 2 , Koi.2l/2 A0O1.21 \ [Koi,2i/2] 2 S c &-cj(A/3oi.2i/2ir) 

TE01-TM2I c H cv K01.21 A0oi.2i 2 K ° I - 2lScwcv(A/3 ° 1 ' 2l/2,r) 



IV. COMPARISON OF APPROXIMATE AND EXACT SOLUTIONS WITH 
MEASURED RESULTS 

In this section the approximate results obtained in the last section for 
TE01 — TM21 mode conversion, the measured results, and the results 
from a numerical integration of the coupled line equations in (7) are 
compared for different curvatures of the waveguide axis. Let us first 
consider an idealized bend, as shown in Fig. 5a, having constant curva- 
ture in the horizontal plane and a sinusoidally varying curvature in the 
vertical plane 

c H (z) = 1/R H 
cy(z) = A.v sin — 2 (19) 

Am 

From Table I, we can approximately model the TM^ conversion 
process with the following set of two coupled equations. Since S CHCV is 
much greater than S C 2 V and S^ H at A/3oi.2i/2ir c/m the TM2T level is much 
greater than the TM 2 i level, which is neglected. 

— Aoi = - jhoiAoi + ;Koi,2i TT~ A v sin ( "T" z ) A 2i 
dz Hh \A m / 

— A 2 \ = jK 0h2 i— A v sin (—z) A 0l - jh 2 iA^ (20) 

dz Rh \A m / 

The coupling term in (20), sin (2ir/\ m )z, is sinuosoidal, and an analytic 
solution for this case has been obtained by Miller. 13 Miller's solution for 
(20) will be compared with both the numerical solution and the sec- 
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Fig. 5 — Comparison of predicted losses for constant curvature bend with sinusoidal 
curvature in orthogonal plane. 



ond-order perturbation theory solution in the following. After the so- 
lutions for Aoi and A^[ in (20) are determined as in Miller, 13 the nor- 
malized TM^ level, G^, follows from (11): 



Gn(L) = 



[ I /Aft \ 21 1/2 1 

C.L sin C. | l + (fj[r) I L 



I /A/3,\2|l/2 



(21) 



where 



C, - -Koi t 2i — Av, 
Z Hh 



2ir 
Aft, = Aftoi,2i - — 

Am 



The corresponding value of Goi is 

|Goi(L)|2 = l-|G 2 -I(L)|2 



(22) 



(23) 



from which the loss is readily obtained. From (14) we obtain the sec- 
ond-order perturbation theory solution for G^ for the two curvature 
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functions in (19) as 

G 2 1(L) = C.L (24) 

In Figure 5b, we compare the results from (21) and (24) with a direct 
numerical integration of the corresponding seven-mode system in (7) 
for the curvature data in Fig. 5a. We see that there is a substantial dis- 
agreement at 100 GHz, the frequency at which A0 O i,2i is equal to 2ir/\ m , 
because the propagation constant of TEoi in a bend is slightly different 
from that in the straight waveguide. This causes the frequency at which 
we obtain coherent coupling to differ slightly from 100 GHz. On per- 
forming the numerical integration at a number of frequencies in the vi- 
cinity of 100 GHz, it was found that a frequency of 101.5 GHz resulted 
in a maximum TM^ conversion as also shown in Fig. 5b. The 101.5 GHz 
curve differs by less than 2 percent from the result predicted by (21). The 
second-order perturbation theory solution in (24) is also in good agree- 
ment for small TM£1 levels, as for C.L < 0.1, C,L « sin C.L. From the 
results in Fig. 5 we draw two important conclusions. The equivalent 
two-mode model in (20) can be summed to give an accurate represen- 
tation of the seven-mode system in (7) and the second-order perturbation 
theory solution (14) is valid for small couplings, C.L < 0.1. 

A controlled experiment was conducted in order to compare the an- 
alytic results discussed above with measured data. An experimental bend 
was constructed having constant curvature in the horizontal plane and 
either "zero or very small" curvature or a sinusoidally varying curvature 
in the vertical plane. The measured curvature for this line is given in Fig. 
6. Shuttle pulse loss 14 measurements were carried out on this bend for 
the two different vertical curvature cases. From the difference of the two 
measurements the loss due to TM^J conversion was determined and is 
shown in Fig. 7 (curve A). Two predicted loss curves are also given in Fig. 
7; curve B from a numerical integration of the coupled line equations (7) 
for the data in Fig. 6, and curve C from the sinusoidal coupling theory 
(21). There is excellent agreement between all three sets of curves with 
the exception of some slight frequency shifts. The frequency shift be- 
tween B and C has been discussed. The shift between A and B is probably 
due to small errors (<1 percent) in our estimated values for the differ- 
ential propagation constant A/?oi,2i* 

A similar comparison has been carried out on a 720 m long section in 
the WT4 field evaluation test. The measured curvature for this section 
is given in Fig. 8. There are substantial changes in profile in this section 
and vertical curvatures greater than 0.020 m -1 were measured in several 
regions. The measured loss characteristic for this section is given in Fig. 
9. A numerical integration of the coupled line equations, for the measured 
curvature data in Fig. 8, was carried out at discrete frequencies, as in- 
dicated by a circle or cross in Fig. 9 for two cases. The first case allowed 
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Fig. 6 — Measured curvature data for experimental route bend. 

only for first order TM]i, TM^, TEl 2 , and TE^ conversion, and it can 
be seen that the predicted losses are much lower than the measured 
losses. The second case allowed for TML TMH, TEl 2f TE^, TMJi, and 
TM2I conversion. It can be seen in Figure 9 that there is excellent 
agreement between the measured loss and the numerical integration 
results. The absolute loss for the peak at 114 GHz is « 1.5 dB. This is a 
fairly large loss, and because of this the perturbation theory will yield 
losses that are too large. The perturbation theory losses are approxi- 
mately 30 percent greater than the numerical integration estimate, as 
also shown in Fig. 9. 

The large measured loss peak at 112 GHz is principally due to TM|i 
conversion occurring at the large vertical curvature peaks at 125 and 250 
m in Fig. 8. This is readily seen in the measured "moving piston" trace 
of Fig. 10. Here we are measuring 11 the reflected TE01 level at the end 
of the line as a perfectly reflecting piston is moved through this mode 
filter section in a left-to-right direction in Fig. 8. The large increase in 
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Fig. 7 — TEoi loss due to mode conversion in modulated bend at Chester. 

the TE i ripple at points A and B corresponds to regions of high vertical 
curvature in Fig. 8. 

The numerical integration used in solving (7) consists of an evaluation 
of the transmission matrix for 5 cm long incremental sections of the 
waveguide line. A predicted moving piston trace is readily obtained on 
multiplying the transmission matrix with its transpose. The final result, 
given in Fig. 11, agrees well with the measured trace in Fig. 10. 

V. LOSS VARIATION IN BENDS 

In the preceding discussion we have seen that TE01-TM21 mode 
conversion effects result in sizable losses for the TE i mode in bends. 
This will have a significant impact on repeater siting and spacing for a 
commercial system and it is desirable to determine the sensitivity of this 
loss component to system parameters such as guide size (a), frequency 



TE 01 TRANSMISSION WITH AXIAL CURVATURE 1863 



0.012 




800 



Fig. 8 — Measured curvature for 720-m-long section in WT4 field evaluation test. 

(/), bend radius (Rb), rms curvature (cvrnJ, etc. On examination of the 
loss components in Table I, we see that the TM21 loss depends on the 
product of a coupling coefficient #01,21 and a geometrical factor S CHCV 
orS c 2_ c 2. The coupling coefficient K"oi,2i is the product of two curva- 
ture coupling coefficients, e.g., Coi.n and Cn >2 i which vary as af, and the 
inverse of a differential propagation constant, e.g., AjSqui. which varies 
as a 2 /. Thus /f 01,21 varies approximately as a 4 / 3 - The TM21 conversion 
loss is proportional to If 01,21 > and varies as a 8 / 6 . 

In a plan bend with an arc length MOO m and a constant radius Rb 
< 100 m, the geometrical factor S CHCV (A@oi t 2i/2ir) is approximately equal 
to 1/Rb 2 S cv (A l Soi i 2i/27r). The magnitude of S CHCV is much greater than 
in nominally straight waveguide. Thus the TE i loss due to TM^ con- 
version is much greater in a plan bend than in straight waveguide. 
Conversely, the value of S c * - c » (A/?oi,2i/2ir) is increased to a lesser degree 
from the value in straight waveguide for spatial frequencies greater than 
0.04 c/m and thus the increase in TE01 loss due to TM^ conversion is 
of less significance. 

With the above assumptions the additional TE01 loss per unit length 
in a plan bend may be approximately written as 
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Fig. 9— Loss comparison for 720-m section. 

^-«*^©W©'(?fei) dB ' m < 25 > 

Here a is the guide radius in mm, / is the frequency in GHz, R B the bend 
radius in m, and S cv the vertical curvature power spectral density in 
l/m 2 /c/m. As an example, from (25) a 95 m radius bend with a net angle 
of 75 degrees or an arc length of 131 m has an average added loss of 0.12 
dB at 110 GHz assuming a value of 7 X 10~ 5 l/m 2 /c/m for S cv . This value 
of S c v is typical of the WT4 field evaluation test. 10 

Figure 12 is the measured curvature for a bend in the WT4 field eval- 
uation test with nominal parameters similar to those in the preceding 
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Fig. 10 — Moving piston measurement of relative TE i level for 720-m-long section at 
113 GHz. 



paragraph. The TM21 mode conversion loss for this bend was calculated 
directly using the result of perturbation theory from Table I in Section 
IV. This agrees fairly well with the numerical integration of the coupled 
line equation as shown in Fig. 13. Several bends from the field evaluation 
test were analyzed and exhibited similar agreement. The perturbation 
theory losses were generally 10 to 20 percent higher than the losses 
predicted from the numerical integration. The fine structure of the loss, 
with the exception of small shifts in the loss peaks, was similar for both 
cases. Figure 13 also gives a comparison of the average added loss pre- 
dicted from (25) and the actual loss for the bend. In general (25) is useful 
only for estimating the average loss for a given bend. The actual loss for 
a given bend exhibits significant fluctuation about this average value 
as shown in Fig. 13. Any practical installation will contain several route 
bends in a given repeater span. These bends are independent and thus 
an expression similar to (25) is useful as it gives an accurate estimate of 
the average added bend loss. 



VI. CONCLUSION 

A coupled line system of equations has been presented which predicts 
the effect of second-order TM21 mode conversion on TE01 loss. A nu- 
merical integration of the coupled line equations developed from this 
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Fig. 11— Calculated TEoi level for 720-m-long section at 113 GHz. 



model agreed very well with measured results for a bend experiment 
designed to have high TM21 conversion and also agreed with measure- 
ments on sections of the WT4 field evaluation test. A second-order per- 
turbation theory solution for the TE01 loss due to TM21 conversion was 
derived and discussed. The perturbation theory results agreed well with 
the results obtained by numerical integration of the coupled line equa- 
tions for typical sets of curvature data from the WT4 field evaluation test. 
The perturbation loss estimates were approximately 20 percent higher 
but the two loss predictions had the same fine structure except for slight 
shifts in the frequencies of the loss peaks. 

An approximate formula for the additional loss due to route bends was 
given. The added loss is proportional to the inverse square of the bend 
radius. The added loss increases as the eighth power of guide radius and 
the sixth power of frequency. The added loss is approximately propor- 
tional to the mean squared vertical curvature. This approximate formula 
works well on the average but does not predict the fine structure of the 
loss for a given route bend. 

There are several areas worthy of further investigation. Additional 
work is required on the statistics of the TM21 mode conversion loss, 
particularly the loss statistics for a single bend or pair of bends. We have 
also assumed, without proof, that the TMu and TE12 loss is given by the 
first-order perturbation theory 4 and is little affected by TM21 conversion. 
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Fig. 12 — Measured curvature for plan bend in WT4 field evaluation test. 

Note that significant energy passes through these modes to the TM21 
mode and an examination of higher-order terms is desirable. Finally, 
the three-mode set of equations (10) has great significance for a practical 
dielectric-lined waveguide system and is of as much interest as the 
two-mode systems investigated to date. Improved approximate solutions 
for systems of this type would be desirable. 
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APPENDIX 

Derivation of coupling coefficients 

In many cases of interest 123 circular waveguide with a complex jacket 
structure may be replaced with a simpler "wall impedance" model with 
a resultant simplification in the analysis. This is true for the WT4 
waveguide and thus we have the following set of equivalent boundary 
conditions at the inner waveguide surface r = a: 

Zq = £ /rY z | r=a 
Z z = — E z /H,p\ r = a 

The total fields within the guide can be expressed as a superposition of 
two scalar functions. 

T p = AV p (x„r) sin p<f> 

T\, = N ll J p (x ll r) cospcf) 

Here v and p are integers and the eigenvalue x v a is defined by 

k " = XA (42) 

The fields vary as e~J hvZ in the axial direction and k v is the solution of 

the characteristic equation 

jwe aZ z 

ksJ p {k v ) J' p (k„) + jweoZ&Jpik,) 

h* r i =0 (43) 

-^~«/|(W + ^p(M[cr p (^)+7«6oZ |fJp(MJ 

The total fields are expressible as a superposition of normal mode 
fields 

E = £ At(z)ete-J h >* + £ A;(z)e;e + J h >* (44) 

V V 

H = E At(z)hie-J h >* + E A;(z)h;e + J h >* 

V V 

where 

uta \ or rd(f>/ coeo \ro<j) or I 

+ z(±)j^ \^^ 2 T v (45) 
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+ 4>{±)VT-(— + d v — — ) +z;o)€oV — d—T v 
The separation constant d v is given by 

a " kiV* kw j 

p h 

y» = J' P (k v )lkJ p (k v ) 

The modes are orthonormal as the vector cross product 



= -— 2 (y» + J/"«>aZ z ) (46) 



J f e+xhf - ds - ± «, M (47) 

The normalization constant N„ is given by 
jJVJMJJC*,) [ (l + ^) (l - ^ + 2y„ + fcjyj) 

for p^O, and by 



*N*k*Ji(k v ) -ff ((1 + fcjyj + 2y„)} = 1 (49) 



djh? 
k- 
for TEo„ modes, and by 

aftgftVKMa + *W + 2y,)} = 1 (50) 

for TMo„ modes. 

For a nonideal waveguide where the inner surface is of the form r = 
a + tpi((f>,z) where e « 1, we may express the total field in the distorted 
guide (E,H) as a superposition of the normal modes 5 in an ideal guide 
of radius a. On taking Maxwell's equations, multiplying by /ij or ej, 
integrating by parts over the guide cross section and adding we find 



C C (VXE=-ju>noH)-htds+ C C (VXH = j«€oE).c*ds 
reduces to 

■f Aj(z) = ±-( (f EXhZ-fidl+ (f HXe^-tdl] (51 

OZ 2 I tSr=a ^ K? r=a 
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We now apply the boundary condition in (1) at r = a + epi(<f),z). For 
any initial fieldEg^Ho in undistorted guide the distortion epi yields a 
distorted field E, H to first order of 



E = E + eEi 
At r = a + cpi we have to first order from (40) 



(52) 



1 d Pl dE 0l p dp±zj . „ dH Qz . w ] r „\ 

_ T7 A 0r+Pl~T^ + A 1*~ Z * ^~"0r+Pl~T — + «lz (53) 



a d<f> "' ■ " dr 



dz dr 



§£} rp gggg „ _ „ I^Plrr ggog „ 

iS r + Pl „ +£,i 2 --Z 2 /ior+Pl . +"1* 

dz or a d0 dr 

From (51) we need to determine 

[EXh; + HXe;]\ r=a -f 

On substituting [from (52)]_for E and H and using boundary condi- 
tions (40) at r = a for Eq and Ho we find 

[EXh; + HX el\ ■r\ a = e[E 1 Xh; + H 1 X ej] • f\ a (54) 

The right-hand side of (54) may be reduced to terms containing only Eq, 
Ho, and /ij by expanding and substituting for {E\^ — Z<pHi z ) and (E\ z 
+ Z Z H U ) from (53). 
On letting 

Eq = EA„e„ and H = £^A7 

i' J* 

we finally obtain a system of coupled equations of the form 

f A J(z) = £ ff£' +) A+Cz)e>'< fc *- ft '>« + K<:-- , A;(2)e^+' J ^ (55) 

Only forward coupling is significant as h v » /i,, and we need concern 
ourselves only with the forward coupling coefficient, K^; +) which is of 
the form 

K - "2 Jo " ^La. e " P+Pl ar + Ma^^ + Pl aJJ 

+ ^["a5 6 -" Pl ^ +Z *(^ + Pl ^)]) fld * (56) 

This is similar to Unger's 15 result for helix waveguide (Z^ = 0) coupling 
coefficients but contains several additional terms. 

K*f is the coupling coefficient for displacement of the waveguide walls 
from a right circular cylinder. If pi(0,z) is of the form 
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epi((f>,z) = e cos 0/(2) 

we can determine the curvature coupling coefficient (C*;, + ) from K*'* 
as in Rowe and Warters. Since the curvature 1/R(z) is approximately 
equal to —ef w (z) we have 

C<:: +) = Kfr +) /\j(h, - h v )f (57) 

This is the normalized curvature coupling coefficient referred to in the 
main text of this paper. 
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